Hadronic dissipative effects on elliptic flow in ultrarelativistic heavy-ion collisions 
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We study the elliptic flow coefficient V2(rj, b) in Au+Au collisions at ^J~s~ = 200 A GeV as a function 
of pseudorapidity n and impact parameter b. Using a hybrid approach which combines early ideal 
fluid dynamical evolution with late hadronic rescattering, we demonstrate strong dissipative effects 
from the hadronic rescattering stage on the elliptic flow. With Glauber model initial conditions, 
hadronic dissipation is shown to be sufficient to fully explain the differences between measured i>2 
values and ideal hydrodynamic predictions. Initial conditions based on the Color Glass Condensate 
model generate larger elliptic flow and seem to require additional dissipation during the early quark- 
gluon plasma stage in order to achieve agreement with experiment. 

PACS numbers: 25.75.-q, 25.75.Nq, 12.38.Mh, 12.38.Qk 



One of the important new discoveries made at the Rel- 
ativistic Heavy Ion Collider (RHIC) is the large ellip- 
tic flow v-2 in non-central Au+Au collisions [l|. At the 
highest RHIC energy of *fs = 200 A GeV, the observed 
V2 values near midrapidity (|t7|<1), for not too large 
impact parameters (b<7 fm) and transverse momenta 
(pT ^1.5 GeV/c), agree with predictions from ideal fluid 
dynamics 0] , including 0, the predicted dependence 
of V2 on the transverse momentum px and hadron rest 
masses JjJ. From these observations it has been con- 
cluded [6j that in these collisions a quark-gluon plasma 
(QGP) is created which thermalizes on a very rapid time 
scale Ttherm < lfm/c and subsequently evolves as an al- 
most ideal fluid with exceptionally low viscosity. 

On the other hand, the ideal fluid dynamical descrip- 
tion gradually breaks down as one studies collisions at 
larger impact parameters and at lower energies Q or 
moves away from midrapidity 0, 0, 0, This has 

been attributed alternatively to incomplete thcrmaliza- 
tion of the QGP during the early stages of the expansion 

ipative 

viscosity"). It has recently been argued LTJJ that 
quantum mechanics imposes a lower limit on the shear 
viscosity of any medium, but that the shear viscosity 
of the QGP can not exceed this lower limit by a large 
factor [18|. On the other hand, qualitative arguments 
were presented in Ref. |l6j which emphasize the impor- 
tance of hadronic dissipation and support a picture of a 
"nearly perfect fluid strongly coupled QGP (sQGP) core 
and highly dissipative hadronic corona" in ultrarelativis- 
tic heavy-ion collisions. The importance of viscous effects 
for a successful description of RHIC data on V2 and i>4 
was also emphasized in 0], although this work left it 
open whether the corresponding lack of thermalization 
occurs mostly at the beginning or towards the end of the 



[13J ("early viscosity") and/or to dissi pati ve effects dur- 
ing the late hadronic expansion stage [lj, lla, 1 161 ( "late 



expansion phase. In the present paper we explore this is- 
sue more quantitatively, by trying to answer the question 
how much of the observed deviation of i>2 from the ideal 
fluid prediction can be attributed to "late viscosity" in 
the dissipative hadronic phase, and whether or not signif- 
icant additional dissipative effects during the early QGP 
stage are required for a quantitative understanding of the 
data. 

Our study is based on a comparison of a hybrid model, 
combining an ideal fluid dynamical QGP stage with a re- 
alistic kinetic description of the hadronic stage (hadron 
cascade), with data on the centrality and rapidity de- 
pendence of the pt- integrated elliptic flow V2{f],b) for 
charged hadrons 1201 . We find that with Glauber model 
initial conditions |2l| . suitably generalized to account 
for the longitudinal structure of the initial fireball [22^ . 
hadronic dissipation is sufficient to explain the data. On 
the other hand, initial conditions based on the Color 
Glass Condensate (CGC) model (2^, 13 ^c&d to larger 
elliptic flows which ovcrpredict the data unless one ad- 
ditionally assumes that the early QGP stage possesses 
significant shear viscosity, too, or that the QGP equa- 
tion of state is significantly softer than usually assumed. 
Our analysis points to a need for a better understand- 
ing of the initial conditions in heavy-ion collisions if one 
hopes to use experimental data to constrain the QGP 
viscosity and equation of state. 

A (l+l)-dimensional hydro+cascade model was first 
proposed in Ref. |2^ . putting emphasis on radial flow 
in heavy-ion collisions. It was later extended to 2+1 
dimensions for the study of elliptic flow near midrapid- 
ity fl4ll26j |. By combining a hydrodynamic description of 
the early expansion stage with a hadron transport model 
at the end we can implement a realistic treatment of the 
freeze-out process and of viscous effects during the fi- 
nal hadronic phase. Here we extend the above models to 
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full (3+l)-dimcnsional hydrodynamics [27j, in order to be 
able to study the rapidity dependence of elliptic flow. Let 
us briefly summarize our model. For the hydrodynamic 
stage, we solve the conservation laws d^T^ = with the 
ideal fluid decomposition T^ v = {e+p)u ,1 u v —pg^" (where 
e and p are energy density and pressure and m m is the fluid 
4-velocity) in Bjorken coordinates (r, x±,rj s ) We 
neglect the finite (but at RH1C energy very small) net 
baryon density. A massless ideal parton gas equation of 
state (EOS) is employed in the QGP phase (T > T c = 1 70 
MeV) while a hadronic resonance gas model is used at 
T <T C . When we use the hydrodynamic code all the way 
to final decoupling, we take into account [inj ] chemical 
frcczcout of the hadron abundances at T c h = 170 MeV, 
separated from thermal freezeout of the momentum spec- 
tra at a lower decoupling temperature Idee, as required 
to reproduce the experimentally measured yields |28| . 

For the hydro+cascade description, a hadronic trans- 
port model JAM \2$k is employed for the late stage of 
the expansion. JAM simulates nuclear collisions by in- 
dividual hadron- hadron collisions. Soft hadron produc- 
tion in hadron-hadron scattering is modeled by exciting 
hadronic resonances and color strings. Color strings de- 
cay into hadrons after their formation time (t~ lfm/c) 
according to the Lund string model PYTHIA j^l • Lead- 
ing hadrons which contain original constituent quarks 
can scatter within their formation time with other 
hadrons assuming additive quark cross sections |3lj . In 
the current study, it is initialized with output from the 
above (3+l)-dimensional hydrodynamics by using the 
Cooper-Frye formalism |32^ (rejecting backward going 
particles) [lj,|26|. We switch from hydrodynamics to the 
cascade approach at the switching temperature T sw = 169 
MeV, i.e. just below the hadronization phase transition. 

We here study two types of initial conditions for the 
evolution. The first, which we call "modified BGK initial 
condition" [22I Y\?\ , assumes an initial entropy distribu- 
tion of massless partons according to 

T^tn-M)/^) 



drj s d 2 xj_ 
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where ?y s = i ln[(i+z)/(t— z)] is the space-time rapidity, 
x±_ = {x, y) is the position transverse to the beam axis, 
C = 24.0 is chosen to reproduce the measured charged 
hadron multiplicity in central collisions at midrapidity 
|34) . Y(, is the beam rapidity, and f pp is a suitable 
paramctrization of the shape of rapidity distribution in 
pp collisions, 
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with A?7=1.3 and 0^ = 2.1, which are so chosen as to 
reproduce the measured j_>seudorapidity distributions for 
charged hadrons j3^ 



N ££ and N " 



arc the number of 



wounded nucleons in the two nuclei and the number of 
binary nucleon-nucleon collisions, respectively, as calcu- 
lated from the Glauber model nuclear thickness function 
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with the inelastic nucleon-nucleon cross section a 
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42 mb and r± = \_{x±^b) 2 + y 2 ] 1 ^ 2 (where b is the im- 
pact parameter). The soft/hard fraction a = 0.85 was 
adjusted to reproduce the measured centrality depen- 
dence |34j of the charged hadron multiplicity at midra- 
pidity. At ?7 S = 0, Eq. 0) reduces to dS/(drj s d 2 x±) cx 
[a(A^ art +A^ art ) + (l-a)AUi]/(l+a) 0; this param- 
eterization is equivalent to the one used in Ref. [23j, 
~ ^(N£ H +N£ lt ) + xN mlh with z=^§. From 
Eq. IjTJl. we can compute the entropy density at the ini- 
tial time to = 0.6 fm/c Q of the hydrodynamic evolution, 
s(tq, x±, rj s ) = dS / (TQdi] s d 2 x±), which provides the ini- 
tial energy density and pressure distributions through the 
tabulated EOS described above. 

The second typ e of initial conditions is based on the 
CGC model |3J. Specifically, we use the Kharzeev- 
Levin-Nardi (KLN]approach [23j in the version previ- 
ously employed in |24j . In this approach, the energy dis- 
tribution of produced gluons with rapidity y is given by 
the fey-factorization formula |38( 
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d 2 px 



d 2 x±dy N 2 - 1 J p T J 4 
x <Pa(xi, {p T +k T f /^,x±) 
x (f>B(x2, (pr-k T ) 2 /4:;x±). 



-a s {Q 2 ) 



(6) 



where x\_2 =pt exp(±y) / 'y/s and pt is the transverse mo- 
mentum of the produced gluons. We choose an upper 
limit of 3GeV/c for the pr integration. For the uninte- 
grated gluon distribution function we use 



b A (x,k 2 ,;x±) = 



kC f Qs 
2ir*a s (Q-i) Q'i+A' 2 ' 
kC f Ql 



k T < Qs 

k T > Q s 



(7) 



where C f : 



n: 



?N — and Q s denotes the saturation momen- 
tum. We introduce a small regulator A = 0.2GeV/c in 
order to have a smooth distribution in the forward ra- 
pidity region \y \ > 4.5 at RHIC (other regions are not af- 
fected by introducing this small regulator). The overall 
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normalization n is determined by fitting the multiplicity 
of charged hadron at midrapidity at ^/snn = 200 GeV 
for the most central collisions. The saturation momen- 
tum Q s of nucleus A in A+B collisions, needed in the 
function <f>A, is obtained by solving the following implicit 
equation at fixed x and x±: 



Q 2 s (x,x ± ) = ^a s (Q 2 s )xG(x,Q 2 s ) 
Of 



d 2 xj_ 



(8) 



An analogous equation holds for the saturation momen- 
tum of nucleus B in For the gluon distribution func- 
tion G inside a nucleon we take the simple ansatz [2^ 



xG(x,Q 2 ) = K\n 




(9) 



with A = Aqcd = 0.2 GeV. We choose K = 0.7 and A = 0.2 
so that the average saturation momentum in the trans- 
verse plane yields (Q^(x=0.01)) ~2.0GeV 2 /c 2 in central 
200 A GeV Au+Au collisions at RHIC. For the running 
coupling constant a s in Eq. JSJ we use the standard per- 
turbative one-loop formula, but introducing a cut-off in 
the infra-red region of small Q s (i.e. near the surface of 
the nuclear overlap region where the produced gluon den- 
sity is low) by limiting the coupling constant to ct s < 0.5. 
We can obtain the energy density distribution at time To 
from Eq. (BJ, e(ro, x±, rj s ) = dEx/(Todr) s d x±), where y 
is identified with r) s , and use this as the initial distribu- 
tion for the hydrodynamic evolution. 



In Fig. ^ we show the initial eccentricity 
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(y2 +x 2) 01 

the source at midrapidity (rj s = 0) for our two models for 
the initial conditions. Here (• • • ) represents the average 
taken with respect to the initial energy density distribu- 
tion e(ro, x±, rj s = 0). While the BGK model interpolates 
between the binary collision and wounded nucleon scaling 
curves (being closer to the latter), CGC initial conditions 
are seen to give much larger eccentricities. This can be 
traced back to a steeper drop of the energy density profile 
near the edge in the CGC model. In ideal fluid dynam- 
ics, the larger eccentricities e result in larger elliptic flow 
coefficients V2 

This is shown by the thin lines in Fig. [3 which compare 
ideal fluid dynamical calculations with RHIC data from 
the PHOBOS collaboration [2c|. Note that the hydrody- 
namic calculations shown here use an EOS in which the 
hadron abundances are held fixed below T c at their chem- 
ical freeze-out values established during the hadroniza- 
tion process For this EOS it is known from earlier 
hydrodynamic studies 0, [Ijj that the slope of pion px 
spectrum stalls during the hadronic evolution. Without 
an additional transverse velocity kick already at the be- 
ginning of the hydrodynamic stage it reproduces the 
pion data from central collisions at midrapidity only up 
to ~1 GeV/c, falling off too steeply at larger px- In con- 
trast, the proton pt spectrum continues to flatten dur- 




FIG. 1: (Color online) Initial spatial eccentricity e= 
at midrapidity as a function of impact parameter b, for 
200 A GeV Au+Au collisions with CGC (solid red) and BGK 
(dashed blue) initial conditions. For comparison we also show 
initial conditions where the initial parton density at midra- 
pidity scales with the transverse density of wounded nucleons 
(dotted green) and of binary collisions (dash-dotted black) 



ing the hadronic stage, due to continued build-up of ra- 
dial flow, reproducing the proton data up to ~2GeV/c 
if one chooses Td ec = 100 MeV as the freeze-out tempera- 
ture. The hybrid hydro+cascade model gives almost the 
same proton px spectrum as ideal fluid dynamics with 
Tdec = 100 MeV, but further flattens the pion spectrum 
at larger px, thereby improving the agreement with the 
pion data up to ~2 GeV/c. This can be understood as a 
consequence of shear viscosity in the hadronic phase, with 
the fast longitudinal expansion generating positive addi- 
tional viscous pressure components in the transverse di- 
rections whose effects on the px spectrum grow quadrat- 
ically with px flij and decrease with increasing hadron 
mass. The agreement with the measured px spectra is 
slightly better for CGC than for BGK initial conditions. 

As explained in the chemically frozen hadronic 
EOS results in smaller 1/2 values than for a hadronic EOS 
which assumes hadronic chemical equilibrium all the way 
down to kinetic freeze-out at Td ec - Correspondingly, our 
curve for BGK initial conditions lies below the one shown 
in Fig. 2 of |§3] which uses the latter EOS. It is notewor- 
thy that for central and semicentral collisions the data 
seem to lie consistently somewhat above the hydrody- 
namic predictions with BGK initial conditions. While at 
first sight this seems to argue against the validity of the 
BGK model, it must be noted that event-by-event fluc- 
tuations in the geometry of the nuclear overlap region, 
which we don't take into account tend to signifi- 
cantly increase the measured V2 in collisions with small 
impact parameters (large iV par t values) |45j . 
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FIG. 2: (Color online) pr-integrated elliptic flow for charged 
hadrons at midrapidity (| r\ j< 1) from 200j4GeV Au+Au 
collisions, as a function of the number JV par t of participating 
nucleons. The thin lines show the prediction from ideal fluid 
dynamics with a freeze-out temperature Td cc = 100 MeV, for 
CGC (solid red) and BGK (dashed blue) initial conditions. 
The thick lines (solid red for CGC and dashed blue for BGK 
initial conditions) show the corresponding results from the hy- 
dro+cascade hybrid model. The data are from the PHOBOS 
Collaboration 



The difference between the eccentricities given by the 
BGK and CGC initial conditions may seem surprising in 
view of the fact that the centrality dependence of hadron 
multiplicities in both approaches can be made numeri- 
cally very similar by a proper choice of parameter a in 
Eq. [23j. When parameterized in terms of Eq. QJ, 
the main prediction of the CGC approach is the near 
independence of a on the collision energy which is con- 
firmed by the data. The reason for the big difference in 
the eccentricities stems from the different entropy pro- 
files predicted by the two approaches, especially in the 
regions where the density of produced particles is rela- 
tively small. While these differences contribute little to 
the total observed multiplicity, they appear quite impor- 
tant for the evaluation of eccentricity. The application 
of the CGC approach in the region of small parton den- 
sity is of course questionable, and a better theoretical 
understanding of the transition from high to low density 
regimes is clearly needed. 

The thick lines show the centrality dependence of 
V2 from the hydro+cascade hybrid model. Whereas 
with BGK initial conditions the dissipative effects of the 
hadronic phase in the cascade model reduce the purely 
hydrodynamic i>2 sufficiently to bring the theory in agree- 
ment with the data even for peripheral collisions, CGC 
initial conditions, with their larger eccentricities, cause 
so much elliptic flow that even the hybrid model over- 
predicts the data significantly. This is true for both the 
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FIG. 3: (Color online) The pseudorapidity dependence of V2 
for charged hadrons in (a) central (3-15%), (b) semicentral 
(15-25%), and (c) peripheral (25-50%) Au+Au collisions at 
y/s = 200 A GeV. The corresponding impact parameters are, 
respectively, 6 = 4.0, 6.3, and 8.5 fm. The hydrodynamic evo- 
lution is initialized with modified BGK initial conditions. The 
lines show the predictions from ideal fluid dynamics with 
T dcc = 100 MeV (solid blue) and T dcc = 169 MeV (dashed 
green). The red circles show the corresponding results from 
the hydro+cascade hybrid model. The black squares are mea- 
surements by the PHOBOS Collaboration l2(il[ . 



px-integrated elliptic flow (shown in Fig. 0) and its pr- 
slope, dv2{pr) / dpr , at midrapidity, for both pions and 
protons (to be published elsewhere). With CGC initial 
conditions even the QGP phase must exhibit significant 
dissipative effects if one want to reproduce the RH1C 
data. 

In Figs-EOandlUwe show the measured [2(| rapidity de- 
pendence of vi for three centrality classes, together with 
calculations at three representative impact parameters 
6 = 4, 6.3 and 8.5 fm corresponding to these centrality se- 
lections (i.e. adjusted to give the correct average number 
of participants A?p art in each case as quoted in Ref. [2fj|). 
Results from BGK initial conditions are shown in Fig. [3] 
The lines show ideal fluid dynamical calculations with ki- 
netic decoupling assumed at T^ ec = 100 MeV (solid blue) 
and Tdcc = T sw = 169 MeV (dashed green). The dashed 
lines undcrpredict the data at all impact parameters and 
all but the most forward rapidities, indicating the need 
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FIG. 4: (Color online) Same as Fig. El except for using CGC 
instead of BGK initial conditions. 



for generating additional elliptic flow during the hadronic 
stage below T c . The solid lines, on the other hand, 
strongly overpredict the forward rapidity data in scmipe- 
riphcral and peripheral collisions, showing that ideal hy- 
drodynamics generates too much additional elliptic flow 
during the hadronic stage. The hydro+cascade hybrid 
model (red circles) gives a good description of the data 
over the entire rapidity range for all three centralitics, 
with the exception of the midrapidity region in the most 
central collision sample as already discussed above. The 
hadronic cascade model provides just the right amount of 
dissipation to bring the ideal fluid prediction down to the 
measured values, especially in very peripheral collisions 
and away from midrapidity. 

The situation is different for the more eccentric CGC 
initial conditions, as shown in Fig. 0] Now v-i is over- 
predicted at all centralities and rapidities if the hadronic 
phase is described by ideal fluid dynamics, and the dissi- 
pative effects of the hadronic cascade arc no longer suf- 
ficient to reduce v-i for the more peripheral bins enough 
to obtain agreement with the data. In the peripheral 
sample (25-50%) the excess elliptic flow persists at al- 
most all rapidities and exists even if the fireball freezes 
out directly at hadronization (no hadronic evolution at 
all). Significant dissipation in the early QGP phase is 
needed in this case to bring the theoretical prediction in 



line with experiment. 

We conclude that the answer to the question, whether 
all of the observed discrepancies between elliptic flow 
data and ideal fluid dynamical simulations can be blamed 
on "late hadronic viscosity" and fully eliminated by em- 
ploying a hydro+cascade hybrid model such as the one 
studied here, depends on presently unknown details of 
the initial state of the matter formed in the heavy-ion 
collision. With BGK initial conditions hadronic dissipa- 
tion seems to be able to reduce the elliptic flow enough 
to bring the theoretical predictions in line with the data, 
leaving little room for additional dissipative effects in the 
early QGP stage. CGC initial conditions yield signifi- 
cantly more eccentric sources and produce larger than 
observed elliptic flow even if dissipative effects in the 
late hadronic stage are taken into account. In this case, 
the early QGP phase must cither be significantly softer 
than parametrized by our equation of state or exhibit 
significant viscosity itself. With CGC initial conditions, 
the excess over the V2 data from Au+Au collisions at 
RHIC persists even at midrapidity in all but the most 
central collisions; with such initial conditions, the stan- 
dard claim that at RHIC energies the measured elliptic 
flow at midrapidity exhausts the theoretical upper limit 
predicted by ideal fluid dynamics must be qualified. 

We see that a data-based attempt to establish lim- 
its on the viscosity of the quark-gluon plasma requires, 
among other things, a better understanding of the ini- 
tial conditions of the fireball created in RHIC collisions. 
Unfortunately, very few direct probes of the initial con- 
ditions are available. In Ref. [2^, E3| 3-dimensional jet 
tomography was proposed to test the longitudinal struc- 
ture of BGK and CGC initial conditions in noncentral 
collisions. A specific feature of CGC initial conditions is 
a predicted sign flip of the first Fourier moment v\ of nu- 
clear modification factor Raa(pt, y,<P) at high px as one 
moves away from midrapidity | 46| . Alternatively, one can 
try to exploit the fact that the large (even if not perfect) 
degree of thcrmalization observed in heavy-ion collisions 
at RHIC limits the amount of entropy produced during 
the expansion. When taking into account that final state 
rescatterings in the medium produce only very small ef- 
fects on the shape of the rapidity distribution, the finally 
observed charged hadron rapidity distributions therefore 
severely constrain the initial entropy and energy density 
profiles (2j|. A better theoretical understanding of the 
initial conditions, especially of the transition from the 
high density to small density regimes, is needed to ex- 
tract the viscosity of quark-gluon plasma at the early 
stages. A systematic study of the charged hadron ra- 
pidity distributions for a variety of collision centralities, 
center of mass energies and system sizes is needed to as- 
sess which description of the initial state yields a more 
consistent and efficient overall description of all available 
data. 
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